Parallel heterodyne detection of dynamic light scattering spectra from gold 
nanoparticles diffusing in viscous fluids. 
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We developed a microscope intended to probe, using a parallel heterodyne receiver, the fluctu- 
ation spectrum of light quasi-elastically scattered by gold nanoparticles diffusing in viscous fluids. 
The cutoff frequencies of the recorded spectra scale up linearly with those expected from single 
scattering formalism in a wide range of dynamic viscosities (1 to 15 times water viscosity at room 
temperature). Our scheme enables ensemble-averaged optical fluctuations measurements over mul- 
tispeckle recordings in low light, at temporal frequencies up to 10 kHz, with a 12 Hz framerate array 
detector. OCIS ; 180.6900, 090.1995, 290.5820 
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Advances in dynamic light scattering microscopy [l]- 
01 have recently raised the possibility of developing pas- 
sive microrheology tools based on the measurement of 
local thermally-driven intensity fluctuations Q . The 
related homodyne measurement schemes involve usu- 
ally single- mode detection tTriEj- Concurrently, hetero- 
dyne detection is a powerful method for investigating 
dynamic light scattering (DLS) in low-light conditions 
[lol \vi\ . In usual heterodyne detection schemes, the 
intensity fluctuations of the mix between the scattered 
field and a local oscillator (LO) field are measured on a 
single-mode detector. Since the probed scattered field 
matches the LO beam mode structure [12], such het- 
erodyne schemes cannot perform multimode measure- 
ments. Multimode homodyne techniques based on in- 
tensity fluctuations measurements have been developed 
to take advantage of the spatial bandwitdh of array de- 
tectors in order to perform speckle ensemble averaging 
and/or spatially-resolved measurements l3-17|. How- 
ever, these parallel homodyne techniques suffer from a 
temporal frequency bandwidth imposed by the camera 
frame rate and from limited sensitivity in low light. 

We present here a DLS microscope designed for prob- 
ing the first-order Doppler spectrum of the light scattered 
by nanoparticles in Brownian motion. Our microscope 
performs heterodyne detection with a camera used as a 
multi pixel detector. The Doppler spectrum is explored 
by sweeping the LO beam optical frequency, so that the 
temporal frequency bandwidth is not limited by the cam- 
era frame rate. 

The optical configuration is sketched in FiglTJ The 
whole optical apparatus has been described elsewhere 
[isj . A laser beam (Al = 532 nm, field i?L, angu- 
lar frequency ojl, average power 50 mW, single axial 
mode, CW) is split into tunable-frequency local oscil- 
lator (LO) and illumination arms (fields _Elo and Ei) to 
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FIG. 1: Optical configuration. A colloidal gold suspension in 
a calibrated viscous fluid is illuminated by single mode laser 
light. The total internal reflection at the top interface ensures 
a dark-field illumination of the particles. The light field E 
scattered by the sample beats against a local oscillator fleld 
-Blo, detuned by Aoi with respect to the illumination fleld 
El. Ei^o is tilted by with respect to E. The interference 
pattern \E + i?Lo|^ is recorded by a CCD camera. 



form a Mach-Zchnder interferometer. The samples are 
thin chambers (thickness h 100 fim, several millime- 
ters wide) containing an aqueous solution of gold parti- 
cles (diameter d = 40 nm : BBI international, EM.GC40) 
with added sucrose to increase the viscosity. The particle 
density is p ~ 2.0 x 10^^ m~^. At Al = 535 nm, the inte- 
gral scattering cross-section cr^ of a 40 nm gold particle 
normalized by the particle volume v is as/v ^ 5.2 ^m.~^ 
[l§|. Since the optical mean free path is I — {pas)~^, 
the ratio of the sample thickness to the mean free path 
is thus about h/l ~ 3.5 x 10"''', optical fluctuations are 
interpreted with single scattering formalism in dilute sus- 
pensions [l3| • A Dove prism guides the illumination field 
El to the sample. Microscope immersion oil between 
the prism and the slide enables refractive index match- 
ing. To prevent the collection of the unscattered field, 
a total internal reflection is provoked at the glass-air in- 
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FIG. 2: (a) Hologram (1024 x 1024 pixels) calculated in the 
exit plane of the objective (arrow), (b) Zoom (141 x 141 
pixels) of the objective pupil. 



terface on top of the sample. The incident field makes 
an angle a ~ 50° in water with respect to the optical 
axis. The field scattered by the gold beads E <C i?LO is 
collected by a microscope objective (Leica PL Fluotar: 
100 X magnification, NA = 0.75, air, exit pupil diam- 
eter ~ 3 mm). The LO linear polarization is tuned to 
match the dominant polarization axis of the object field. 
Off-axis optical mixing (A6' = 0.08 rad tilt) of E and 
E\,o with a beam splitter (BS) results in an hologram 
/ — \E + Ei^of' that is recorded with a CCD camera 
(PCO Pixelfly QE, 1392 x 1024 pixels of d 6.45 ^m, 
frame rate ws/(27r) = 12 Hz, exposure time te — 83 ms) 
from which only the central N^xNy = 1024 x 1024 pixels 
region is used for the calculations. 

The DLS spectra are acquired by sweeping the detun- 
ing frequency Aw of the LO. At each frequency point 
Aw, 16 successive CCD frames {/i, ...,/i6} are recorded. 
The field component E in the detector plane is calculated 
from the difference of consecutive images In+i — In- This 
quantity is back-propagated numerically up to the exit 
plane of the microscope objective by a standard convolu- 
tion method [23| with a reconstruction distance D = 34.9 
cm. A representation of the squared amplitude of the 
back-propagated hologram is shown in Figl5](a). Because 
of the angular tilt Ad between the interfering beams, the 
image of the objective is shifted in the upper right corner 
of the hologram. Figl2] (b) shows a zoom of the objective 
pupil image. The heterodyne signal si oc \E\^ is averaged 
within the square region shown in Figl2] (48 x 48 pixels 
with pixel size XD/{Nxd) — 28 ^m). 

In the time domain, the normalized autocorrelation 
function oj^ of the fluctuating object field in the detector 
plane is [10| : 



9id{q,T) = exp( 



-Dq^r) = exp(-r/Tc), 



(1) 



where Tc is the relaxation time, q is the scattering vector 
(momentum transfer), and D is the particle diffusivity. 
Since the beads undergo Brownian motion, the average 
mean-square displacement of each particle at time t scales 
as ~ Dt, where the diffusion constant D is linked to the 
thermal energy /cbT and to the mobility /i of the parti- 



cles: D = fiksT. In the limit of low Reynolds number 
fi — l/(67r ryr) , where ry is the viscosity of the medium 
and r is the (hydrodynamic) radius of the spherical parti- 
cles. The illumination and collection geometry introduce 
momentum transfer diversity, and because we seek to 
study dynamic fluctuations in non-varying optical scat- 
tering conditions, we retain 



gidir) « exp{~Dak'^T) 



(2) 



where k — 27r/AL is the angular wavenumber and a is 
an adjustment parameter used to describe the scattering 
structure factor of the optical configuration. Addition- 
ally, the electric field E that reaches the detector is a 
superposition of a field Eg statically scattered and a field 
Ed dynamically scattered by the sample 2l|, \2^. The 
normalized autocorrelation function gi of the total field 
E = Ed + Es is 



5i(t) = (1 - p) l5id(T)| +p, 



(3) 



where p = \Es\ /{\Ed\ + \Es\ ). According to the 
Wiener-Khinchin theorem, the first-order power spec- 
tral density si(w) of the field fluctuations is the Fourier 
transform of the field autocorrelation function gii^r). Let 
sid{u}) be the Fourier transform of |(7id(T)|. Taking into 
account the instrumental response B in the frequency 
domain, we have : 



si{uj) (1 - p) sirf(w) * B{uj) + pB{uj), 



(4) 



where * is the convolution product. In the absence of dy- 
namic scattering, the off-axis component of the recorded 
image /„ at time tn = 2n7r/a;s is the interference term 
EElQexp{iAujtn)- The measured signal 7°' 
grated over an exposure time te is : 
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(5) 



where i'^ — —1. The normalized apparatus lineshape 
B{Alu) = llfr,^'''' - /,f -"'^•1V(16|7;P|£;loP) is then : 



B{Aoj) = (AwTE)"^sin2(AwTE/2) srn^ {ttAlj/ojs}- (6) 

This function is plotted on fig. [3] (line 5). Hence, in 
good approximation, frequency shifts Aw probe actual 
fluctuations frequencies w. The normalized first-order 
power spectrum of the fluctuating field is the Fourier 
transform of equation [T] : 



Sid(w) = l/(l+wV2), 



(7) 



where Tc = Girrjr / (q^ksT). Experimentally, the samples 
must be prepared with great care to avoid stray light. 
Otherwise, the estimation of the cutoff frequency t^^ is 
made difficult because both fluctuating and non fluctuat- 
ing light are recorded within the available CCD dynamic 
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are obtained from a preparation without Brownian parti- 
cles; they are in agreement with the apparatus lineshape 
derived from EqlHl 

These results demonstrate the suitability of the pro- 
posed method for parallel heterodyne dynamic light scat- 
tering measurements of nanoparticles dynamics in micro- 
scopic volumes. The resolution, range and sensitivity of 
first-order spectrum measurements potentially enable a 
robust exploration of a wide range of viscosity-dependent 
fiuctuations at microscopic scales. 
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FIG. 3: First-order light fluctuations spectra normalized by 
the relative weight of dynamically scattered light. The rela- 
tive viscosities of the samples [23] are r;/?7o = 1-00 (1), 1.69 
(2), 3.94 (3), and 15.40 (4). Continuous lines are plotted 
from EqlH Data points (5) represent the apparatus lineshape 
in the frequency domain, and the continuous line is plotted 
from Eq|^ The horizontal axis is the detuning frequency 
Aa;/(27r). 

range. A high relative weight of static light prevents an 
accurate detection at high frequencies because it reduces 
the effective dynamic range. 

The measured spectra from samples of relative viscos- 
ity ?;/77o = 1.00 (1), 1.69 (2), 3.94 (3), and 15.40 (4) are 
represented in FiglSl where rjo = 10~^ Pa.s is pure wa- 
ter viscosity. For each measurement, the frequency is 
detuned from 3 Hz to 10* Hz over 80 points with log- 
arithmic spacing. Each scan lasts ~ 106 s and 5 min- 
utes are needed to save the data. The frequency scans 
are performed twice with 15 minutes intervals and sig- 
nal variability from one measurement to another is visi- 
ble at 10^'^ from the maximum typically (points clouds). 
The si lines plotted in Figl3] are calculated from Eq|4j 
with ksT = 4.11 X 10"^-^ J at room temperature (298 
K), r = 40/2 = 20 nm, and ^ 3jfc2^ The adjust- 
ment parameter a — 3.7 has been chosen to get a good 
agreement with the experimental data points. The ex- 
pected momentum transfer from particle scatterin g of 
the incident and reflected beams, in DLS conditions [10| 
is = 4n^ sin^(^/2) fc^. = 1.33 is the refrac- 

tive index of water. For the incident beam, we have 
= 1.3 fc^ {i3 ^ a = 50°). For the reflected beam, 
we have q"^ = 5.8P {l3 = n - a = 180 - 50°). The de- 
rived cutoff frequencies are Tc^^ — 888 Hz (1), 524 Hz (2), 
225 Hz (3), 58 Hz (4); these cutoffs are inversely propor- 
tional to the viscosity and spread over a wide frequency 
range. The values of the relative weight of the dynamic 
light component for the different preparations used to fit 
si{u!) to the experimental data are 1 — p = 0.48 (1), 0.15 
(2), 0.14 (3), 0.13 (4), 0.0 (5). Data points (5) of FigE] 
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